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Abstract 

The construction of all 2-D Lagrangian systems, which admits besides 
the energy another integral of motion that is quartic in velocities was 
reduced in our previous article (J. Phys. A: Math. Gen., 39, 5807-5824, 
2006) to a single nonlinear PDE. 

In the present note we introduce a new solution of this equation, lead- 
ing to a new integrable system with a quartic integral, which involves 16 
free parameters. A special case of the new system admits interpretation 
in a problem of rigid body dynamics. It gives a new integrable variation 
of the cases due to Kovalevskaya (1889), Chaplygin (1903) Goriatchev 
(1916) and Yehia (2006). 



1 Introduction 

1.1 Historical 

The famous Kovalevskaya's integrable case of rigid body dynamics was the first 
example of a mechanical system that admits an integral of motion quartic in 
velocities [1 . For more than a century this case attracted attention of many 
specialists, who treated explicit solution in terms of time and gave several mod- 
ifications and generalizations. Only in the last three decades there appeared a 
few new integrable systems with a quartic integral, but mainly concerning the 
motion of a particle in the Euclidean plane under the action of certain potential 
forces. A short, but nearly complete up to its date, list of those systems can 
be found in Hietarinta's review paper [2 . A few cases of the same type were 
obtained later in [3J. 

In virtue of Maupertuis principle, the motion of a natural mechanical system 
can be brought into equivalence (in the orbital sense) with the geodesic flow on 
some Riemannian metric. Metrics on the Riemannian sphere associated with 
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known integrable cases of rigid body dynamics were constructed in [3]. Two 
families of integrable systems with a quartic integral on S 2 were obtained in [5] 
and [6]. Few more works discussed possible integrable systems with low-degree 
polynomials on S 2 and the hyperbolic plane H 2 (see e.g. [7]-[H]). 

The method introduced in our work |10) and developed in several later works, 
has led to construction of a large number of several-parameter families of in- 
tegrable systems with a complementary integral ranging from second to fourth 
degree (see e.g. [T3], [IS])- This method leads in a natural way to several- 
parameter integrable systems that allow as special cases, for particular values of 
the parameters, interpretations as motion on various flat and curved 2-D man- 
ifolds. Although our primary interest is in systems on Riemannian manifolds, 
some of the constructed integrable systems live on pseudo-Euclidean or pseudo- 
Riemannian manifolds. New examples of integrable rigid body dynamics are 
common by-products of this method [Tt)]-[T7]. 

The culmination of this method was the construction of the so-called " mas- 
ter" system with a quartic integral [16] . It involves the largest ever number of 21 
parameters and covers almost all systems of that type that were known earlier. 
Applied to rigid body dynamics, this system gave rise to 4 new integrable cases 
of motion under action of potential forces. In two cases the potential is a single 
valued function on the configuration space, but multivalued in the other two 



The problem of constructing all integrable mechanical systems admitting a 
quartic complementary integral on a Riemannian manifold has been reduced 
in [16 to a single nonlinear PDE, which we call the resolving equation. A 
solution of this equation is what we call generating function. This generates 
a conservative lagrangian system together with its quartic integral, valid on 
its zero-energy level. In the present article we introduce a new solution of 
the resolving equation. We also formulate a theorem that has been implicitly 
used in our previous work to construct systems integrable on arbitrary energy 
levels from condtional ones usually built on the zero-energy level. In its final 
unconditional form the integrable system based on the new solution involves 
15 parameters. Special cases of it introduce new integrable cases of the types 
investigated earlier by Toda, Bozis and others for a particle moving in the 
Euclidean plane. A special case of the new system adds a rare new integrable 
problem in rigid body dynamics. 

1.2 Construction of integrable systems 

According to a result of Birkhoff [18], the general natural mechanical system 
(on an arbitrary 2D Riemannian configuration space) can always be reduced in 
certain (isometric) coordinates £, r/ and time parametrization r to the form of a 
ficticious plane system described by the Lagrangian 
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so that the equations of motion become 

x „ = dU dU 
dx dy 

are restricted to their zero-energy level 

£' 2 + r/ 2 - 2U = (2) 

The energy constant h for the original system enters linearly as a parameter in 
the ficticious force function U, which has the structure 

U = A(h~V) (3) 

where V is the potential of the original system and A is a function that depends 
on the metric of the configation space. 

According to a result of [TD] (see also [H]), if an integral of motion of the 
mechanical system exists in the form of a polynomial of the fourth degree in ve- 
locities, this integral can be reduced while preserving the form of the Lagrangian 
([T]) to the form 

I = C' 4 + PC' 2 + QfY + R = Io (const) (4) 
where R is given by the quadrature 



dU 

" = - / Q-az*n- 



in which []q means that the expression in the bracket is computed for r\ taking 
an arbitrary constant value 770 (say), and the other three functions involved are 
expressed as 

in terms of a single generating function F which satisfies the nonlinear partial 
differential equation 

d 2 F d A F d 4 F d 3 F d 3 F d 3 F d 3 F 
d 2 F d 4 F d 2 F d A F 

= (7) 

The set of solutions of this equation generates all systems of the type ([1]) having 
an integral of the form (Q} on the zero level of their energy integral. Affecting 
all possible conformal mappings of the complex £ = £ + ir\ plane followed by a 
general point transformation to the generalized coordinates q± , qi with a suitable 
change of the time variable we obtain all systems of the general form on 2D 
Riemannian (or pseudo-Riemannian) manifolds, having a quartic integral on 
the zero level of their energy integral. 
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1.3 Certain forms of solutions of the resolving equation 



The fourth-order PDE ([7} is homogeneous nonlinear of the second degree. It is 
not known whether this equation is integrable, in the sense that some procedure 
can be pointed out to construct all its solutions, nor it can be reduced to an 
equation of lower order. The best one can do is to try to construct solutions 
involving as much arbitrary parameters as possible. 

1. Equation (J7J admits a solution of the form 

F = A*(0*(»?) (8) 
where A is a constant and both $ and \& satisfy the same equation 

* ( %) + 5 ^ ^ - 6o4*(0 - (9) 

in which 04 is a separation constant. Comparing with the results of |16) 
we conclude that the original isometric variables £, r\ are not practically 
suitable coordinates for the description of the solution, and that a sym- 
metric separation solution (J5J can be more conveniently expressed in the 
form 

F = Apq 

in terms of the pair of variables p, q related to £, rj by the relations 



dp 



\J a,4p 4 + a 2 p 2 + aip + a 

v = r . dq (10) 

J y/atq* + b 2 q 2 + b x q + 6 

The three constants 03,01,00 (and 62,61,60) are integration constants of 
the fourth-order equation © for $ (and The fourth constant is elim- 
inated by a shift of the variable p (q) . 

A quite rich system, first announced in [16] as the master system, can be 
derived from the solution 



in which 

Fo = 

Gq — 



F = F (p) + G (q) + vpq (11) 

dp f p (4C + 4Cip + 4 Ap 2 + \b 3P 3 )dp 



{/a^ + a 3 p 3 + a 2 p 2 + dip + a J (aiP 4 + a 3 p 3 + a 2 p 2 + a x p + a ) 3 / 4 
dq f q (4D + AD iq + AAq 2 + \a 3 q 3 )dq 



t/aiq 4 + b 3 q 3 + b 2 q 2 + hq + b J Kg 4 + hq 3 + b 2 q 2 + b x q + 6 ) 3/4 

(12) 
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and the relations 

dp 



\J a 4 p A + a 3 p 3 + a 2 p 2 + a\j> + a 



v = r , dq as) 

J ^Ja 4 q A + b 3 q 3 + b 2 q 2 + hq + b Q 

generalize (ITU)) by a shift of the two variables p, q. This system has the 
Lagrangian 



1 r P 2 q 2 



L 



2 yj a 4 p 4 + a 3 p 3 + a 2 p 2 + (Hp + a \J a 4 q A + b 3 q 3 + b 2 q 2 + b\q + b 
^ vb 3 p 3 + Ap 2 + dp + C + ;/a 3 g 3 + Aq 2 + D iq + D . 
y a4P 4 + a 3 p 3 + a 2 p 2 + a\p + a y/ a±q A + b 3 q 3 + b 2 q 2 + hq + bo 
q{Aa A p 3 + 3a 3 p 2 + 2a 2 p + ax) p{4a 4 q 3 + 3b 3 q 2 + 2b 2 q + b 



\J a 4 p A + a-3P 3 + a-iP 2 + aiP + a-o \J a 4 q 4 + b 3 q 3 + b 2 q 2 + hq + bo 

which involves 15 free parameters ao,ai,a 2 ,a 3 , 04, bo, b%, b 2 , b 3 , A, Co, Ci, Do, D\, v 
and the quartic integral valid only on the zero level of its energy. This in- 
tegral was provided in [Jj5] . It turned out that the integral can be reduced 
to somewhat simpler form 

{p 2 /2 + vq{Aa 4 p 3 + 'ia 3 p 2 + 2a 2 p + ai ) + vb 3 p 3 + dp + Co} 2 . .. 

I = j 3 5 wpq 

a^p + a 3 p a + a 2 p z + a\p + a 

-%v 2 \J 'a 4 p i + a 3 p 3 + a 2 p 2 + a x p + a \J a 4 q A + b 3 q 3 + b 2 q 2 + b~iq + b 

-Av{v[b 2 p 2 + a 2 q 2 + 3pq(2a 4 pq + b 3 p + a 3 q)} + 2Apq + D lP + dq}5) 



2 Real dynamics: 

2.1 Classes of time-parametrized geodesies. 

To construct a system that would be a priori integrable on all levels of energy 
we should receive a function U from © that has the structure (J3J), i.e. involving 
in a linear way an arbitrary parameter h. In that case The generating function 
should have the form F = hFo + F 4 , where Fq,F\ are functions not depending 
on h. This is usually expected when solving linear PDEs. However, we shall 
rely on solutions of the nonlinear PDE (JJJ having this property of being a linear 
superposition of solutions. It is obvious from that the function F can always 
be reduced to the forms XF (A arbitrary constant), but in practice we seek F 
in the form of a linear superposition of such expressions, so that 

fc 

i=i 
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where Xi,i = 1, ...,k are arbitrary parameters and {Fi} do not depend on any 
one of the A$s. This choice will lead to 

U = J2x i U i ,U i = -\(^ + ^) (17) 

and {f/j} do not depend on any one of the AjS. Each Ai can be regarded as the 
energy parameter of the real mechanical system before reduction to the form (fTJ) , 
and hence we call them energy-like parameters. However, for a much general 
result we formulate the following 
Theorem 1 

Let the system with the (reduced) Lagrangian 

k 

in which ()' = -4- and {Ui} do not depend on any one of the Xj, admit on its 
zero- energy level an integral 

i = HZ,v,t',v,*i,-~,*k) (is) 

then the Lagrangian 



E frUi 



L * = \(Z «<^)K 2 + v 2 } + '-f (20) 

1=1 E oaUi 

i=i 

in which () = 4r and t is the natural time, admits an unconditional integral 
( valid for arbitrary initial conditions ) of the form 

r=I(£,ri,Ai,Mi,aih + 0i,-..,a k h + p k ) (21) 

k 

where a±, a k , Pi, Pk are arbitrary parameters, A = E ct%Ui and h is the 

i=i 

energy integral for i20\) . The energy parameter h can be either substituted 
by its value taken on an arbitrary motion of the svstem ![W)\) or replaced by its 
functional form in the state space 



h 



E Pi u i 

i=l 



i=l 



i=i 



Proof 
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Introduce new arbitrary parameters {a;,/?,} by the substitution = ath + 
Pi, i = 1, k. The reduced force function in (| 18[) splits into the two expressions 

k k 

^2 PiUi + h ^2 OLiUi. Now, we make the change of the independent variable r 

i=l i=l 

to the original (natural) time variable according to the relation 

k 

dr=^-,K = Y j a i U i (22) 

i=l 

It is easy to see that those substitutions lead to the Lagrangian 

Li = V + h 

The additive parameter h in Li is insignificant and discarding it reduces Li to 
L*. The zero level of the energy integral for Li is the same as the /i-level for L* . 
Notes: 

1. Although we are concerned here with systems with a quartic integral, in 
this theorem the integral (fl~9f is a general function of its arguments and 
not necessarily quartic nor even a polynomial in velocities. 

2. According to Maupertuis principle (see e.g. [20]), it should be noted that 
the whole 2k— parameter family of mechanical systems (|20j) on their level h 
of energy share one and the same fc— parameter geodesic flow (or the same 
fc— parameter family of unparametrized trajectories) on the manifold with 
metric 



fc 

ds 2 = 2(^2 \ l U l )[d£ 2 + dr] 2 } (23) 



3. The form (1181) in the independent variable r can be regarded as the 
simplest form of the 2fc— parameter family of mechanical systems. Let 
(£(t),T)(t)) be the general solution of (fT5|) . The general solution of (f2"0"j) 
can be expressed in terms of the natural time t by first integrating the 
relation (|22|) to get 

t= f Adr= fj2 a i u Mr),v(r))dr (24) 

J i=l 

and solving for r we obtain r = r(t), so that the solution of ([2"U|) is 
(£(t(£)), /y(r(t))). The form of the function r(t) for each choice of the set 
of parameters {a^} determines the way in which the mechanical system 
(|20j) . which lives on the -generally speaking- non-Euclidean manifold with 
metric 

k 

ds 2 = 2C^a l U i )[di 2 +dr) 2 ] 

i=i 

describes its trajectory in the plane of the isometric variables £, r\. 
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3 A new solution of the resolving equation 



In the present section we try a deformation of the master system to accomodate 
on more term of the product type, probably at the expense of enforcing certain 
restrictions on some of the parameters figuring in (|12[) . In fact, we assume F in 
the form 

F = F + K lP q + K 2 p 2 q 2 (25) 

Inserting the expression (|2"5"|) in equation ([7]) and using the relations (|13[) . we 
obtain after some manipulations a polynomial expression in the two variables 
p, q that must vanish identically. Those equations have been solved under the 
condition that 04 7^ 0, so that the two polynomials that occur under the fourth 
degree root signs are of the fourth degree. This resulted in only one new case, 
which differs from the master system of [16 . This case will be considered in 
detail in the next subsections. 



3.1 The generic restricted case 



We write the Lagrangian and the comlementary integral for this case after some 
transformation to more symmetric form that does not affect the generality of 
the system (arbitrary shifts of the two variables will retain the raw case) : 

L=-[ - ^ + - ^ — }+TT (26) 

2 y'au 4 + kiu 2 + ko \ / av i + m\v 2 + too 

_ N[(4av 2 + TOi)(4au 4 + 3kiu 2 + 2k )} + vuv(2au 2 + k x ) + Ku 2 + D 



2y/au 4 + k x u 2 + k Q 
N[(4au 2 + k x )(Aav 4 + 3m x v 2 + 2to )] + vuv(2av 2 + mi) + Kv 2 + E 



(27) 



I = 



2v 'av 4 + m\v 2 + too 

{u 2 + N(4av 2 + mi)(4au 4 + 3km 2 + 2k ) + vuv(2au 2 + k x ) + Ku 2 + D} 2 

{au 4 + k 1 u 2 + k ) 



-4(8Nauv + v)uv - 2(8Nauv + v) 2 \Jau 4 + kiu 2 + k \/ av 4 + miv 2 + m 
-32aN 2 [u 2 (2au 2 + fci)(6au 4 + 3mi« 2 + m ) + k v 2 (2av 2 + mi)] 
-16Na{Eu 2 + Dv 2 + 2Ku 2 v 2 ) 
—4:Nvuv{8a(miu 2 + k\v 2 + 3au 2 v 2 ) + 3fcirni] 

— ^ 2 (toiU 2 + k x v 2 + 6au 2 v 2 ) - AvKuv (28) 

in which the prime represents derivative with respect to the independent variable 
r (ficticious time). Note that when the two parameters ./V and v vanish the 
Lagrangian (|26|) degenerates into a separable one and the integral (f28|) into the 
square of a quadratic integral. 
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The system with the Lagrangian (|26l) admits the integral (|2"51) only on the 
zero-energy level of this system 

1 [ ^ 2 + ^ 2 1 - U — 

2 v 'sji 4 + fciu 2 + fco V 'av A + m\v 2 + mo 

This system involves 10 parameters a, fco, fci, mo, mi, -D, E, K, v 1 N , of which the 
first five enter in both the kinetic energy and potential terms of the Lagrangian 
and the last five ones enter only in the potential terms and, moreover, they enter 
only linearly. The last four parameters constitute a set of energy-like parameters 
(see e.g. [E]). They are essential in building the unrestricted integrable system. 
Comparing ([2"B)) to its counterpart in the "master" system [TB], we find that 
(|26[) involves only one new parameter N , which is not present in the master 
system. When N is set equal to zero (|2^|) turns out to be a special case of the 
master system resulting from one restriction on the coefficients of each of the 
two fourth-degree polynomials entering under the root sign and two restrictions 
on the energy -like parameters (namely, the vanishing of the parameters C\ , D\ 

ofnsi). 

3.2 Dynamics- The unrestricted generalization 

We now proceed to use those parameters to construct a general integrable sys- 
tem valid on arbitrary energy level out of the restricted one. Introducing new 
parameters by the relations 

D = h\ — ha\, E = h,2 — hai-, 
K = hs — ha3, is = hi — ha^, 

N = h 5 -ha 5 (29) 

and performing the change of independent variable to the actual-time parametriza- 
tion by using the relation 

dt 

dr=- (30) 

where 



a\ + a^u 2 + a4uv(2au 2 + ki) + a^miti 4 + v 2 (4au 4 + 3kiu 2 + 2fco)] 

y/au 4 + k\u 2 -I- ko 
j + &3V 2 + a4uv(2av 2 + mi) + a^[kiv + u 2 (4av + 7>m,\v 2 + 2m^J\^ 

V 'aw 4 + m\v 2 + mo 

we arrive at the new Lagrangian 

1 it 2 v 2 

I = -A [ - — + ] — V + h (32) 

2 \f au A + k\u 2 + fco \/ av 4 + miv 2 + mo 
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_ 1 hi + h 3 u 2 + h 4 uv(2au 2 + k x ) + /i 5 [miu 4 + v 2 {Aau 4 + 3kiu 2 + 2k )} 
A y/au 4 + kiu 2 + k 

^ h2 + h 3 v 2 + hiuv{2av 2 + mi) + h 5 [kiv 4 + u 2 {Aav i + 3miv 2 + 2m^ 

y/av 4 + miv 2 + m Q 

which admits on an arbitrary energy level h the integral resulting from (|28[) by 
the substitutions (|29|l and (|30p. i.e. u — > Au, u — > Aw. The integral will depend 
on the parameters occuring in the Lagrangian and also on the energy constant 
h. The last constant may be substituted by its expression in terms of the coor- 
dinates and velocities to get the final form free of the energy restriction. The re- 
sulting system depends on 16 parameters a, ko, k±, mo, mi, ai, a%, 04, as, hi, hi, h%, hi, /15 
and h, of which the first nine enter in both the kinetic energy (the metric of 
the configuration space) and potential terms of the Lagrangian but the last five 
ones enter only in the potential. 

3.3 Special cases 

3.3.1 Generalization of the cases of Bozis and Wojciechowski 

Let a = ko = tuq = l,ki = mi = —2. Under the coordinate transformation 
p = siny, q = sin 2; the Lagrangian (|32|) takes the form 



1 S3 

-fa + sin x siny + 7(2 cos 2 x cos 2 y — cos 2 x — cos 2 y) H = 1 5— ](x 2 + y ) 

2 cos 2 x cos 2 y 

a + b sin x sin y + c(2 cos 2 x cos 2 y — cos 2 x — cos 2 y) H * 1 — 

« V " & ' cos z x cos 11 ?y 



a + sin a; sin y + 7(2 cos 2 a; cos 2 y - cos 2 a; - cos 2 y) + j^-j + 

When — ^ = 61=62 = we have, after ignoring an insignificant additive 
constant 

L = l:(x 2 +y 2 ) — lb sin x siny + c(2 cos 2 x cos 2 y — cos 2 x — cos 2 y)-| — ^ 1 — -1—1 

2 cos" 1 a; cos 2 y 

(35) 

This system is new. It describes plane motion of a particle in a 4-parameter 
potential. The complementary integral of this system can be written as 

. 2 2c? 1 . , o 2g?2 % /7 >. 

I = [x H - — )(y H — ) — 2cosxcosy(o + 2csmxsmy)a;y 

cos 2 x cos y 

+ cos 2 x cos 2 y(6 -I- 2c sin 2 siny) 2 + 4c(di cos 2 y + g?2 cos 2 2) (36) 

When c = this case reduces to a special version of that found by Bozis [2"T] 
and when c = 6 = the system becomes separable and the integral degenerates 
into the product of two quadratic integrals. A slight variation of the parameters 
in (|34[) to be ki = mi = 2 changes trigonometric functions to hyperbolic (or 



(34) 
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exponential) functions, and thus giving a new system like the type of [53]. The 
analog of (|35|) gives a particle in the potential 



V = b sinh x sinh y + c(2 cosh 2 x cosh 2 y — cosh 2 x — cosh 2 y) H \ 1 ^ — 

cosh x cosh y 

(37) 

with the corresponding integral 

I = (x 2 H g — )(?) 2 H 1 — ) + 2 coshx cosh?/(— b + 2c sinh £ sinh y)±y 

cosh x cosh y 

+ cosh 2 x cosh 2 y(6 — 2c sinh x sinh y) 2 + 4c(di cosh 2 y + di cosh 2 a;) (38) 
In a similar way, one can obtain a mix of the two types by taking k\ = —mi = 2. 

3.3.2 Systems of the Toda type 

If in (|32|) we set a = 1, ko = k\ = mo = m\ = 0, the Lagrangian takes the form 



L = ^(a + ae' 2x + (3e~ 2y + 1 e x+y + 5e^ x+ ^)(x 2 + f) 
h + aer 2x + be- 2 y + ce x +y + de 2( - x+y *> 

(39) 



a + ae- 2x + f3e- 2 v + -fe x +y + 5e 2i ~ x +y) 

and the integral may be written, after using the energy integral to eliminate h, 

as 

I = \ 4 x 2 y 2 + 2A 2 [be- 2y x 2 + ae- 2x y 2 + (ce x+v + de 2x+2y )xy] 

+e 2x+2y (c + de x+v ) 2 + 2d(be 2x + ae 2y ) + 4abe- 2x - 2y (40) 

where A = a + ae~ 2x + (3e~ 2y + -fe x+y + 5e 2i > x+y \ 
A special case of the Toda type is 

L = i(i 2 + y 2 ) - (aer 2x + be- 2y + ce x+y + de 2(x+v) ) (41) 

In two more special cases the configuration space degenerates into a plane. Their 
lagrangians can be written in polar coordinates as done in [16] : 

1/-2 2^ ,4a B Ce 20 + De- 2 \ 
L i = -(r 2 +r 2 9 2 )-{Ar 2 + ^ + , ) (42) 



2 

L 2 = -{f 2 +r 2 e 2 )-i + — + —) (43) 



The last three cases of motion in the plane seem to be new, but their potentials 
are not periodic in 6. 
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3.4 Application to rigid body dynamics 



We now consider the general problem of motion of a rigid body about a fixed 
point under the action of a combination of conservative axisymmetric potential 
forces. The equations of motion for this problem can be written in the Euler- 
Poisson form: 





(C-B)qr 


Bq + 


(A - C)pr 


Cr + 


(B - A)pq 


7i + 973 - rj 2 


= 0, 72 + rji 



dV 


dV 


073 


072 


dV 


dV 


= 73^— 
071 


073 


dV 


BV 


0I2 


Oil 



- P73 = 0, 73 + P72 - 97i = 0, (44) 

where A, B, C are the principal moments of inertia, p, q, r are the components 
of the angular velocity of the body and 71, 72, 73 are the components of the unit 
vector 7 fixed in space in the direction of the axis of symmetry of the force 
fields applied to the body, all being referred to the principal axes of inertia at 
the fixed point. 

The system (|44l) admits three integrals: 

h = \{Ap 2 + Bq 2 + Cr 2 ) + V (45) 

h = A Pll + Bq l2 + Cr l3 (46) 

h - 7i 2 +72 2 +7 3 2 = l (47) 

Equations (|44[) admit an equivalent representation in the Lagrangian form 
(see e.g. |16|). which we write here for a dynamically symmetric body, for which 
B = A. As generalized coordinates we use the Eulerian angles: ip the angle of 
precession around the axis of symmetry of the field, 9 the angle of nutation and 
ip the angle of proper rotation (about the axis of symmetry of the body). The 
components of the vector 7 can be expressed as 

71 = sin 8 sin ip, 72 = sin 9 cos <p } 73 = cos 9 

After ignoring the cyclic variable ip in the sense of Routh on the zero level of 
the cyclic integral I2 = 0, the Routhian of this mechanical system expressed in 
the other two angles 9 and if has the form 

Comparing the structure of this Routhian function to that of the Lagrangian 
(|32[) and recalling the procedure followed in a similar situation in |16) . we get 
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convinced that they become identical only in the case of Kovalevskaya configu- 
ration A = B = 2C. 

In fact, setting a = l,k\ — l,fco — 0, mi = 2, mo = 1, a± = — a§ = 

2 a 

0,a2 = Qi3, and affecting the substitution u = g , v = cos <p and renaming 
the remaining parameters, we get the Lagrangian 



L 
V 



hni , sin 2 



2C[a sin 9 sin i/3 + 6 sin 2 cos(2c9) 



A 



cos z 9 sin 



1 + sin 2 



' cos 2 ip 



(49) 
(50) 



and the integral 
sin 6 9p 4 



I = 



(1 + sin 2 9) 4 
sin 2 9ip 2 



(1 + sin 2 6) 2 
29 p> sin 6 cos 6* cos <p 
(1 + sin 2 J) 

-r^ — + 

sm w cos 4 cc> 



2a sin sin <p + 6(2 — 3 cos 2 



2A sin 2 6 



4/x 



cos 2 6> cos 2 iy5 

(a sin 6> + 26 cos 2 6> sin w) + 9 2 \b cos 2p> H 1 

cos^ tp 



4A 



cos 2 9 cos 2 93 sin cos 2 tp 



4a sin 09 „, . . 1 — 3 sin 2 6.. 

1 26(2 sin 2 9 + — ^ — )] 



sin 2 9 cos 2 ip 



(ia 2 - 2a6sin0sinc?)(l - 2 sin 2 6* cos 2 cp) + ^6 2 sin 2 6»[sin 2 6»(cos4os - 1) +4] 



-2A6tan cos 2c? 



(51) 



To facilitate comparison with other results, we now express the last inte- 
grable case of rigid body dynamics in the Euler-Poisson variables, i.e. as a 
solution of the system (|44| , in the next 



Theorem 2: 

For a rigid body with moments of inertia satisfying A = B — 2C and for the 
potential 

V = 2C[ ail + 6( 7l 2 - 7 2 ) + A + 5 2 —^] (52) 

73 72 

equations are integrable on the level I2 = 0. TTie complementary integral 

has the form 

t ,11 , , 1 Kll-lD-,1 . / 2A7i72, 2 
J = (p -g -071 + 673 5 +(2^0-072 5—) 

7 3 73 

+ A[ 2 (p 2 + a 2 ) 73 2 - 2a 7 i7 3 2 - 2A 7l 2 + 26 + %] (53) 
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This case is new. For comparison we provide a table of presently known 
integrable potentials related to the type (|52|) , which admit a quartic integral 
under the condition A = B = 2C: 



Author- year 
Kovalevskaya [I] 1889 
Chaplygin [23] 1903 
Goriatchev[35! 1916 

Yehia [IS] 2006 



The potential (1521 involves a collection of parameters: a of the Kovalevskaya 
type, b of the Chaplygin type, the Goriatchev parameter A and the parameter 
S figuring in both new cases announced in our work [16] . but the combination 
f5"2"j) is new. 
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